Simulation of a two-dimensional model for colloids in a uniaxial electric field 
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We perform Monte Carlo simulations of a simplified two-dimensional model for colloidal hard 
spheres in an external uniaxial AC electric field. Experimentally, the external field induces dipole 
moments in the colloidal particles, which in turn form chains. We therefore approximate the system 
as composed of well formed chains of dipolar hard spheres of a uniform length. The dipolar inter- 
action between colloidal spheres gives rise to an effective interaction between the chains, which we 
treat as disks in a plane, that includes a short range attraction and long range repulsion. Hence, the 
system favors finite clustering over bulk phase separation and indeed we observe at low temperature 
and density that the system does form a cluster phase. As density increases, percolation is accom- 
panied by a pressure anomaly. The percolated phase, despite being composed of connected, locally 
crystalline domains, does not bear the typical signatures of a hexatic phase. At very low densities, 
we find no indication of a "void phase" with a cellular structure seen recently in experiments. 



I. INTRODUCTION 

Colloidal suspensions, small particles dispersed within 
a second medium, are common in everyday life with ex- 
amples ranging from toothpaste to paint, from milk to 
quicksand [lH^]. An interesting class of suspension is 
the electrorheological (ER) fluid 0-Qj]], where the liquid 
medium and the colloidal particles have different dielec- 
tric constants. An externally applied electric field polar- 
izes the particles, and hence the system can be thought 
of as a system of dipolar particles, the moments of which 
are aligned along the field axis. 

ER fluids and their phase behavior have been stud- 
ied experimentally fl2Hl7ll. theoretically 0, |1] and us- 
ing computer simulation |18l422j |. In Ref [18j . Hynnincn 
and Dijkstra used MC simulation to study the phase 
behavior of hard and soft dipolar spheres over a broad 
range of packing fraction and external field strength. The 
morphologies they saw match those seen experimentally. 
What is observed both in simulation and experiment is 
the formation of chains of particles that in turn form 
sheets and body centered tetragonal (bet) clusters, which 
have been interpreted at low temperature T to be the 
crystal phase coexisting with a gas. 

Also very interesting are the experimental reports of 
a "void phase" , where at very low particle concentra- 
tions and sufficient external field strength a cell struc- 
ture appears with relatively colloid rich walls forming 
the border between colloid poor voids. The structure 
has been reported both in the granular (l6j and Brown- 
ian [T3] regimes, where it appears to be a stable structure. 
A similar but transient structure appears as a result of 
a Rayleigh- Taylor-like hydrodynamic instability during 
sedimentation [23]. 

Both the formation of bet crystallites and the void 
phase occur at higher external field strengths, in the 
regime where particles form long chains along the field 
axis. Our goal is to gain a simpler understanding of the 
structures formed in the system, and so we present a sim- 
plified model in which the system is represented by well 



formed (rigid) chains of dipolar hard spheres of uniform 
length. In this case, the system is essentially two dimen- 
sional, with disks representing chains as viewed down the 
field axis. While this is perhaps not physically realistic, 
it should capture the essential features of the real system 
at high field strength, at least in terms of packing and 
dipolar effects. 

Our desire to use this simplified model stems from 
several factors. First, we believe that the finite extent 
of experimental systems plays a significant role in the 
phase behavior seen. Second, the interesting morpholo- 
gies observed experimentally are seen at high fields, in the 
regime where the system is largely composed of chains. 
It would be convenient to think of these morphologies 
arising from interacting chains of dipolar hard spheres, 
but are these morphologies in fact recovered from inter- 
acting chains? If not, what additional ingredients must 
be present to recover them? 

We believe that a set of chains of uniform length 
provide a reasonable starting point for a coarse-grained 
description of the experimental systems, in which this 
length seems to be between 50 and 100 or so spheres. We 
note that simply using hard rods with a uniform dipole 
moment density would not result in attraction between 
chains at short distances. Once the physics of this sim- 
ple system is understood, progressive detail can be added, 
such as the effect of chain length, chain length polydis- 
persity (to better simulate a sediment), fluctuations in 
chain shape etc. 

The main feature that emerges when reducing the sys- 
tem to one of chains, is that the effective interaction be- 
tween chains as a function of their separation, while de- 
pending on chain length, is generally repulsive at long 
range and attractive at short range. This immediately 
suggests the possibility of clusters that are finite in size 
that reduce the likelihood of bulk phase separation. 

The model embeds the finite extent of the system along 
the field axis within the chain-chain interactions. It thus 
provides a convenient way of approaching the bulk three 
dimensional limit without progressively increasing the 
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number of particles in the simulation. In this work, as 
a starting point and for simplicity, we study a system of 
chains each composed of 50 particles, motived for exam- 
ple by Ref. pj|. 

This paper is organized as follows. In Section II, we 
present the details of our two dimensional model and the 
resulting effective potential. We also discuss the energet- 
ics of clustering and provide details on our Monte Carlo 
simulations, as well as on the descriptive quantities we 
calculate. In Section III we present our results, and in 
Section IV we discuss them, especially with regard to 
phase diagrams pertinent to our system. Finally, we give 
our Conclusions and thoughts for future work. 



II. MODEL AND SIMULATIONS 




FIG. 1: (a) The interaction between two particles of diameter 
a and parallel dipole moments, Pi and P2, separated by a 
distance \r\. 8 is the angle formed between the dipole moment 
and the separation vector r. The external electric field E BX t 
is parallel to the z— axis. |(b)| The preferable configuration of 
the dipoles under a relatively strong electric field. 



A. Interacting chains of dipolar particles 

Our approach to modeling a system of colloidal parti- 
cles in a uniaxial electric field begins with the assumption 
that the field induces point dipole moments within other- 
wise hard spheres of diameter a. As the dipole moments 
are always aligned, the effective pair interaction becomes, 
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is the energy scale of the dipolar interaction, a — (ep — 
es)/(ep + 2es) is the polarization, £5 and ep are the 
dielectric constants of the solution and colloids, respec- 
tively, \Eioc\ = \E e xt \ + \E d ip\ is the magnitude of the 
local electric field with a dominant contribution from the 
applied external field and a negligible contribution from 
the field from the other induced dipoles. As shown in 



Fig. 1(a) 9 is the angle formed between the displacement 



vector between two particles and the field axis, which we 
choose to be along the z aixs. 

At sufficiently small reduced temperature T* — 
kpT/Ap), the dipolar spheres tend to align and form 
chains along the field axis as illustrated in Fig. 1(b) 



Thus, we will model the system as one comprising chains, 
each composed of L dipolar spheres. Colloids within 
a chain are assumed to be fixed with respect to each 
other, while each chain as a whole may interact with 
other chains. 

Within this approximation, the crystal phase available 
to the chains is the body centred tetragonal (bet) struc- 
ture, within which chains are either in a "stacked" config- 
uration, a s in F ig. 2(a)| or in a "staggered" configuration, 
as in Fig. 2(b) While the staggered configuration is en- 
ergetically preferred, and indeed is attractive at small 
chain axis to axis separation d, basic geometrical frustra- 
tion prevents all configurations from being staggered. 




FIG. 2: Two chains interacting at closest distance through 
(a) the stacked interaction and (b) the staggered interaction. 



For the stacked geometry, the potential energy as a 
function of d between two chains is given through Eq. [1] 
by summing all interactions between pairs of colloids in 
different chains, 
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As shown in Fig. 3(a) the interaction is always repulsive, 



increases at small d as L increases, but seems to become 
independent of L for larger d. 

For the staggered geometry, the interaction energy is 
given by, 
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As shown in Fig. |3(b)| for small d the staggered interac- 
tion is attractive, increasing in an approximately linear 
way with L. While the range of attraction also grows 
with L, the attractive part is short ranged with a max- 
imum in U g (d) at d w 1.75a for L = 100. At larger d, 
the interaction is repulsive and quickly becomes indistin- 
guishable from the stacked arrangement. 
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FIG. 4: (Color online) The potential energy of an isolated 
cluster of size mxlasa function of I. 



While this analysis is limited because it neglects inter- 
action with other clusters (i.e., density effects), as well 
as entropic considerations, it does foreshadow the im- 
portance of long, string-like structures as well as more 
compact clusters of finite size. 



Monte Carlo simulations 



FIG. 3: (Color online) Chain-chain interaction energy as a 



function of distance (d) for different chain lengths (L), for 
the stacked dipolar interaction, exhibiting weak repulsion at 
large d, and strong repulsion at small d; and (b) the staggered 
dipolar interaction, exhibiting weak repulsion at large d, and 
strong attraction at small d. The inset in each figure shows a 
close-up of the interaction at short distances. 



One motivation for modeling the system in the way 
we do is to highlight the short range attraction and long 
range repulsion between the basic constituents of the sys- 
tem in the regime where chains are well formed. 



B. Energetics of clustering 

The main impact of the short range attraction and 
long range repulsion is the emergence of energetically fa- 
vored finite size clusters. This phenomenon has been the 
topic of numerous studies on colloidal systems, including 
studies pertaining to equilibrium thermal gels. 

To illustrate the effect for our system, we compose rect- 
angular bet clusters to chains wide and / chains long and 
plot for each to the potential energy of each cluster as 
a function of I in Fig. |4j All chains in the cluster have 
length L — 50 and are in contact with nearest neighbors, 
i.e. d/a = \[Z/2 sa 0.866. As the figure illustrates, the 
potential energy for narrow clusters (m = 1, 2, 3) keeps 
decreasing with length, while for wider clusters (to > 4) , 
there appears an energetically favored length. For square 
to x to clusters (not shown) potential energy increases be- 
yond m = 7. 



The system of interacting chains is essentially two di- 
mensional, and as such, we represent each chain of col- 
loidal spheres as a disk in a plane. To capture the short 
range structure, we label disks with either a 1 or -1. If the 
product of the labels is 1, then the interaction is given by 
U c (d;L), the always repulsive stacked potential. Other- 
wise, the disks interact through U g (d; L), which is deeply 
attractive at short range and repulsive for distances be- 
yond 1.64cr. We choose all chains to be composed of 
L = 50 particles, and we simulate N p = 2500 disks. 

We follow the usual Metropolis MC algorithm to gen- 
erate new configurations at constant area [24|. To model 
the slight chain motion in the z direction that can change 
local packing (staggered and stacked), we choose between 
two MC moves, each with 50% probability: a random 
displacement of a disk, or a random displacement of a 
disk accompanied by a change in label of that disk. The 
simulations are done for a broad range of T* and area 
fraction = N p na 2 / '(4Ls 2 ), where Lb is the square 
simulation box length, and we report quantities for state 
points where the mean square displacement of the disks 
after an initial equilibration period is at least a 2 . We 
truncate the potentials at r c = Lb/2. For the largest 
A = 70%, r c = 26.48cr. 

To charcterize our system, we calculate the radial 
distribution function g(r), the structure factor S(q) = 

(PqP-q), where p$ = Np 1/2 J2iLi ex P ' ^»)> and tlic 
potential energy of the system U(T,A^,), which includes 
the standard tail correction assuming g(r > r c ) 1. 
We aslo calculate the distribution of clusters N(n), i.e., 
the number of clusters of size n. We do not consider 
the degree of local crystallinity in defining a cluster, 
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but rather simply define two disks to be in the same 
cluster if they are within a distance of r n = 1.233er 
of each other [U g (r n ) = 01. We do so in the spirit of 
Toledano and coworkers [25| , who studied clustering in a 
system with competing short range attraction and long 
range repulsion in three dimensions, and used the ap- 
pearance of a local maximum in N(n) for n > to de- 
fine a cluster phase. This maximum is equivalently a 
minimum in the free energy of forming an n-sized clus- 
ter, AF(n) = -\n[N(n)/N(l)}. We also keep track of 
whether or not there is a system-spanning cluster and de- 
fine the percolation line by identifying state points that 
have a 50% chance of containing a spanning cluster. 
We report the pressure in terms of g(r) via, 
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where p is the number density per unit area, g c (f) and 
g g (r) are the partial radial distribution functions for 
stacked and staggered interactions, respectively. We have 
accounted for hard core repulsion, omitted the (density 
dependent) impulsive correction at r c , and have added 
the long range tail correction. In this way we approxi- 
mate the behavior of the untruncated system. 

As mentioned in the Introduction, Ref. [l8j indicates 
that at low T and over a range in density, the crystal 
phase coexists with a rare fluid. As another check for 
coexistence within our system in addition to tracking P, 
we wish to track the reduced number density pj of the 
fluid (non-crystal- like) portion of the system as varies. 

To determine the number N c of feci-like particles, i.e., 
disks in a locally crystalline environment, we emp loy the 
method developed by Frenkel and co-workers [26j, based 
on the local bond-order analysis that was introduced in 
Ref. [23] . For each particle we define the quantity 
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where the sum is carried out over all neighbouring parti- 
cles Nb(i) that are within a distance 1.05ct, Y/ m (r ij ) are 
the spherical harmonics calculated for the normalized di- 
rection vector fij between the neighbours. The unit vec- 
tor fjj is determined by the polar and azimuthal angels 
8ij and <pij, although here in two dimensions 6ij = it/ 2. 
Since our crystallites posses square symmetry, we use 



I = 4. For each pair of neighboring particles i and j, 
we calculate the correlation Y2m=-A ^rn{i)qX m {j) , where 

<Z4m = ?4m(«)/El=-4 l<?4m («) | 2 ] 1/2 and q* is the com- 
plex conjugate of q. If the correlation between two neigh- 
bouring particles is greater than 0.05 or less than —0.05, 
then particles i and j are considered to be connected. If 
a particle has at least 3 connected neighbours, then it is 
considered a crystal-like particle. 

Additionally, since we observe crystallites connected 
through a single bridging particle, we also define a parti- 
cle to be crystal-like if it has only two neighbors that are 
both connected to it and the angle between the separation 
vectors connecting the particle with the two neighbours 
is less than a. The choice of a is somewhat arbitrary, 
and so we report results for a — 103° and a = 115°. 
This procedure labels as crystal-like corner particles of 
crystallites as well as compact squares of four particles. 

To determine pf, we first define an area occupied by 
crystallites A c by associating a square of length a to each 
crystal-like disk and accounting for overlaps. The area 
occupied by the fluid is Af = L B — A c and Pf = {N p — 
N c )a 2 /Af. pf should be constant with respect to A^ 
along an isotherm in a region of crystal-fluid coexistence. 

Finally, to test for the existence of a hexatic phase, we 
report on the oricntational correlation function 



G 4 (r) 



(7) 
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k=l 



where Ojk is the angle formed by the bond between par- 
ticle j and neighbor k and an arbitrary fixed axis. In a 
crystal, Gi(r) tends to a constant as r increases. In a 
hexatic phase, it decays as a power law while in a liquid 
it decays exponentially. 



III. RESULTS 

In Figs. [5] and |51 we show sample configurations for 
different T taken from isochores at A^ = 1% — 70%. At 
A^ = 10%, we see that single disks initially cluster into 
short chains (clusters one particle wide) as T decreases, 
but then form more compact clusters at the lowest T. We 
see similar behavior at A<p = 40%, although the higher 
density encourages longer, more intertwined chains that 
thicken as T decreases. At the lowest T, the system com- 
prises rectangular clusters connected with chains. 

The right columns of Figs. [5] and |5] show the progres- 
sion of the system at low T as A^ increases from indi- 
vidual crystallites, to crystallites connected by chains, 
to crystallites with very few chain-like structures. By 
A r f, = 0.70, there are very few non-crystal-like particles. 

What is clear from these pictures is that cluster-cluster 
repulsion, arising from the cumulative long-range disk- 
disk repulsions, prevents the formation of a bulk-like 
crystalline region. The local crystal-like ordering within 



5 



(a) 



lb) 



(d) 



(h) % I Jf 



(c) 






# 




8 


* 


4 





FIG. 5: Sample configurations for A$ equal to (a-c) 1%, (d-f) 
10%, (g-i) 20%, and (j-1) 30%. For As = 1%, configurations 
are shown for T* equal to (a) 4.0, (b) 1.8, and 



A<f, — 10%, 20%, and 30%, configurations are s 
equal to |(d][g][j)| 4.0, [^itTOl 2.3, and^jM 



(c)]0.4. For 
rown for T* 
1.5. 



FIG. 6: Sample configurations for A^ equal to (a-c) 40%, (d-f) 
50%, (g-i) 60%, and (j-1) 70%. Configurations are shown for 
T* equal to Ka^d^gfj)] 4.0, [(b^e^h^k)! 2.3, andf^cjf - 
1.8. 



clusters induced by lowering T does not result in the 
appearance of a crystal phase. It would be difficult to 
classify the system as being one of phase separated crys- 
tal and fluid, despite expecting the system to be fully 
crystalline at higher densities. 

To categorize the structures seen, we plot the phase 
diagram for the system in Fig. 7(a)| where we have in- 
dicated the percolation line (where we define the perco- 
lating cluster by considering bonded particles regardless 
of local order) and a line where AF(n) first exhibits a 
minimum at n > 0, the "cluster line" . 

The AF(n) are shown in Figs. 8(a)|8(c) for area frac- 
tions well within and near the percolation threshold. For 
— 10%, a minimum in AF(n) for small but finite 
n begins to develop at T* w 1.8, moving to larger n 
as T decreases. The transition from the high T fluid, 
characterized by stringy clusters, to the low T cluster 
phase, characterized by more compact, squarish clus- 
ters, appears to be continuous with no free energy bar- 
rier between monomers and larger cluster at low T. For 
A t jj = 30%, the picture is similar, except that the slope 
of AF(n) is smaller in the cluster regime, and that there 
appears to be a small barrier on the order of fcgT sepa- 



rating monomers from larger clusters. While this would 
lend a first-order flavor to the transition, the minimum 
in AF(n) is local and AF(n) increases monotonically 
beyond it. For A<p — 40%, which crosses the percolation 
threshold as T decreases, no local minimum in AF(n) 
forms for small n. Rather, accompanying percolation, a 
local minimum develops at n on the order of the system 
size. 

The potential energy as a function of T for the iso- 
chores is plotted in Fig. [9] For A^ = 1%, where we ap- 
proach the clustered ground state, the potential energy is 
well described by a two-state model, with the exception 
at low T where the vibrational contribution of 3/2fcsT 
is not accounted for in the model. For A^ = 10% and 
20% our lowest T data appear to be just at or below the 
inflection point in U(T), while for A^ > 30%, we have 
not yet located the peak in the heat capacity. Thus, we 
do appear to have a heat capacity peak accompanying 
the change demarcated by the appearance of an n > 
minimum in AF(n), at least at low A^. 

Pressure isotherms, as obtained through Eq.[5]are plot- 
ted in Fig. [TUJ At packing fractions near the percolation 
transition, the pressure shows an inflection, indicating a 
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FIG. 7: (C olor online) Phase diagrams for compar able sys- 
tems. | (a) | Our system of chains of dipolar spheres, (b) The 
three dimensional system of polarizable spheres in an external 
field studied by Hynninen et al. We have redrawn their 
figure after converting packing fraction to area fraction and 
dipolar strength to reduced temperature. |(c)| Phase diagram 
as adapted from 25] for a colloidal system with short-range 
depletion attraction and long-range electrostatic repulsion. 



compressibility maximum, that grows with decreasing T. 
It is not surprising that extrema in response functions 
appear together with percolation (28j . 

Fig. [TT] shows the packing fraction of the non-locally- 
crystal-like, i.e. liquid-like, portion of the system. As dis- 
cussed in the previous section, there is some ambiguity 
in discerning liquid-like particles from crystal-like parti- 
cles, particularly for particles connecting crystallites. As 
such, we plot curves for two values of a (an angle that de- 
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termines whether particles bridging crystal-like domains 
are crystal- like themselves), and find only a quantitative 
difference in the curves. The pf isotherms become flatter 
as T decreases, which is consistent with an interpretation 
of our cluster line as (a remnant of) the phase boundary 
between the fluid and the solid-fluid coexistence phase. 
At T well below this line, pf(A^) is quite flat. 

In Fig.[H we plot S(q) for A = 10%, 40% and 70%. 
In all cases, ordering on the scale of nearest neighbors is 
closely followed by clustering as evidenced by the peak 
near qa = 1. The location of the cluster peak appears 
not to change much with A^, indicating that the distance 
between cluster centers does vary a great deal with den- 
sity. The high value of the cluster peak at A^ = 10% 




FIG. 9: (Color online) Potential energy per particle as a func- 
tion of T* for A$ equal to (a) 1% (red curve indicates two- 
state model fit), (b) 10%, (c) 20%, (d) 30%, (e) 40%, (f) 50%, 
(g) 60%, (h) 70%. 




FIG. 10: (Color online) Pressure isotherms. The inflection 
coincides with percolation. 



may at first seem to indicate a cluster crystal since in 
two dimensions a peak height greater than four is char- 
acteristic of a crystal [Hj]. However, to draw any such 
conclusion, it would be necessary to divide 5(g) by the 
number of particles in a cluster, which is approximately 
four, yielding an effective peak height of two. Visually, 
we do not see ordering of the crystallites. On the other 
hand, for = 70%, 5(g) reaches a height of four at q 
corresponding to nearest neighbors, nominally signifying 
crystallization. 

In Fig. [0 we plot g(r) for A<f, = 10%, 40% and 70%, 
illustrating the T evolution below, near and above perco- 
lation. At 10%, as T decreases, we see the emergence of 
neighbor peaks at short range, while a long range peak 
emerges near 8c only once the cluster line is crossed. 
For A^ = 70%, near T* = 2.0 and below, the neigh- 




FIG. 11: (Color online) Reduced number density of the liquid- 
like portion of the system (p/) as a function of A$ for different 
isotherms at a = 103° (solid lines) and a — 115° (dashed 
lines). 



bor peaks extend to longer distances. The peak heights 
appear to decrease as a power law, as highlighted by a 
power law curve with exponent —1.8 drawn as a guide 
to the eye in Fig. I13f c). Similar, but less pronounced 
behavior is evident at A$ — 60% and 50% (not shown). 
Fig- H5T b) shows A^ — 40% exhibiting cross-over behav- 
ior with a cluster peak emerging near r = 6erat low T. In 
Fig. QU we show G 4 (r) for A^ = 70% and T = 1.8, the 
state point showing the longest range translational order 
through g{r). However, the orientational order quantified 
by G4 (r) appears to decay exponentially, as shown by the 
curve j4exp [— r/£], with A = 0.928 and £ = 11.67<t, over- 
laid on the data. 



IV. DISCUSSION 

This study was motivated by the desire to better un- 
derstand the phase behavior of a system of polarizable 
colloids under the influence of an alternating, external 
field. The system has been hitherto studied as a bulk 
system of dipolar hard sphere particles with all dipole 
moments aligned along one axis. With our reduction of 
the system to one of chains of finite length, we make ex- 
plicit the emergent short range attraction and long range 
repulsion between the constituents of what can be effec- 
tively thought of as a two dimensional system. This type 
of interaction generically promotes finite clustering and 
suppresses the formation of a bulk phase that would oth- 
erwise arise in a fluid-solid coexistence. 

The repulsive part of our potential becomes weaker 
with respect to the short range attraction as chain length 
increases, disappearing in the infinite chain length limit. 
One can readily confirm that the interaction between two 
parallel lines of uniform and aligned polarization decays 
to zero as the length of the chains increases. In this limit 
of long chains, the only source of interaction (attraction) 
is the discrete nature of the spherical dipolar spheres con- 
stituting the chains. Thus, one can see why for a large, 
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FIG. 12: (Color online) The structure factor S(q) for A$ equal 
to 10%, [(b)] 40%, M\ 70%. 



three dimensional system, bulk fluid-crystal phase coex- 
istence emerges. 

For our system at densities below percolation, the high 
T phase is characterized by a distribution of clusters 
sizes, where monomers are the most common state and 
clusters tend to be string-like. As T decreases, the dis- 
tribution shifts so that the most common cluster size is 
finite, and clusters tend to be more compact. The clus- 
ters themselves do not order into a lattice at the T stud- 
ied. At very low density, we expect the transition to the 
cluster phase to continue to occur at progressively lower 
T as the energetics of clustering competes with an en- 
tropically favored monomeric state. At densities above 
percolation, the homogeneous state coalesces into a sys- 
tem of connected crystallites, which increase in size with 
increasing density. 
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FIG. 13: (C olor onli ne) T he pair c orre lation function g(r) for 
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where the dashed 
1 for T* = 1.8 and 



Snapshots of configurations are qualitatively quite sim- 
ilar to experimental ones [Hj], and our results confirm 
that even for uncharged colloids, dipolar repulsions be- 
tween chains stabilize clusters and coarsening (to the 
bulk) should not be expected beyond some finite equi- 
librium size. 

Having suppressed bulk phase separation, we obtain a 
qualitatively different phase behavior from that obtained 
for the three dimensional study by Hynninen (l8| and 
coworkers. For comparison, we present an adaptation of 
their phase diagram in Fig. 7(b) [ after converting their 
packing fraction to area fraction [Ajj = 3/27y), where r\ is 
the packing fraction. Our T scale is rather different, ow- 
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FIG. 14: (Color online) d(r) as a function of r at A^, = 70% 
and T = 1.8. Red curve shows an exponential fit j4exp [—r/^], 
with A = 0.928 and £ = 11.67cr. 



ing to our precondition that our system comprises fully 
formed and unchanging chains. Despite the lack of full 
phase separation, remnants of first order character ap- 
pear in our cluster phase, namely, a possibly non-zero 
free energy barrier separating the monomeric state and 
stable clusters of a larger size, as well as the flatness of 
the pf(A^) curve. As mentioned in the Introduction, it 
would be interesting to study the approach of our two 
dimensional model to the three dimensional case by pro- 
gressively increasing L. 

Our phase behavior appears to be closer in spirit to 
that of a three dimensional system of particles interacting 
through isotropic competing short range attraction and 
long range repulsion studied by Toledano et dl (25|. In 
Fig. 7(c) we show their phase diagram, again rescaling 



packing fraction. One point to emphasize for our system 
is the importance of percolation, for which we see an 
accompanying inflection in the pressure. 

At densities above percolation, we see what appears to 
be the remnant of bulk crystallization in the formation of 
connected crystallites below a T similar to the T at which 
the (unpercolated) cluster phase appears at lower density. 
While a large nearest neighbor peak in S(q) indicates 
crystallinity, the orientational and translational order do 
not extend in space appreciably beyond the character- 
istic size of crystallites. Interestingly, the peaks in g(r) 
decay with a power law dependence, though with a larger 
exponent (~ 1.8) than what is typical for orientational 
correlation function decay in a hexatic phase [29]. 

Even though the crystallites in the percolated regime 
do not seem to align to form a structure with long range 
order, and this may be understood in terms of the de- 
generacy of the long range energy between two clusters 
with respect to precise orientation, at sufficiently high 
density the system must become fully crystalline. Pre- 



sumably, the size of the crystallites continues to increase 
with density, perhaps forming an ordered phase dotted 
with voids, but the nature of the transition is for us an 
open question. 

The transition to the crystal would perhaps be better 
rigorously studied through another model. One draw- 
back of the present model is that it is difficult to handle 
the long-range nature of the interaction in a precise way, 
and the dynamics become quite slow at high density. Our 
potential, which we determine numerically, converges too 



slowly to 



for us to be able to take advantage of 



Ewald techniques. Perhaps the formation of the crystal 
phase could be studied with simpler potentials, or even 
potentials with finite repulsive tails, as were studied re- 
cently in Ref. [3(J. Such a potential would also simplify 
a more thorough study of the thermodynamic behavior 
near percolation. 

Finally, we find no evidence at very low density of the 
emergence of the cellular structure of the "void phase" of 
Ref. [lTj • Clearly, the physics of dipolar hard spheres as- 
sembled into chains monodispersed in length is not solely 
responsible for the appearance of this phase. Perhaps the 
formation of chains helps stabilize structures that begin 
to form as a result of hydrodynamic forces, but this is 
mere speculation. 



V. CONCLUSIONS 

In this work, we reduce a model of aligned dipolar hard 
spheres to one of interacting finite chains of particles, 
treating the chains as a collection of disks in two dimen- 
sions. The long range repulsion between disks (chains) 
promotes the appearance of a cluster phase at the ex- 
pense of bulk crystal formation. Percolation at higher 
density is accompanied by a pressure anomaly. The per- 
colated phase at low T is characterized by connected crys- 
tallites exhibiting orientational disorder that increase in 
size with increasing density and power-law decay of g(r). 
The description of the transition to the crystal phase at 
densities higher than those presented here is of great in- 
terest to us, but such a study would be better done with 
a more computationally tractable potential. We find no 
evidence for the appearance of the "void" phase found 
experimentally at very low density. 
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